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■ A scalar potential coupled to other fields of large disparate masses will exhibit power suppression 

of the quantum loop corrections from these massive fields. Quintessence fields in the dark energy 
regime and infiaton fields during inflation often have a very large background field value. Thus any 
other field with its mass dependent on the quintessence/inflaton background field value through a 
^s^j ' moderate coupling will become very massive during the dark energy /inflation phase and its quantum 

corrections to the scalar effective potential will be suppressed. This concept is developed in this 
paper using the decoupling theorem. The problem then reduces to a quantitative question of the 
' size of suppression effects within the parameter space of coupling constants, scalar field background 

value and renormalization scale. Some numerical examples are presented both for inflation and 
quintessence, but the approach is general and can be applied to any scalar field effective potential. 
The consequences to dark energy of the decoupling effect developed here is that the quintessence 
field need not just be an incredibly weakly interacting field, often added as simply an add-on to 
generate dark energy and have no other purpose. Instead, this quintessence field could play a central 
role in the particle physics dynamics at early times and then simply decouple at late times before 
the onset of the dark energy phase. For inflation a consequence is coupling of the infiaton to other 
fields can be much larger in certain models, without needing supersymmetry to control quantum 
corrections. 
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I. INTRODUCTION 

Often in cosmology, scalar fields with appropriate potentials and interactions are studied in 
order to explain different physical effects. Very well known examples are inflation models, which 
provide a solution for the horizon and flatness problems of the standard cosmology, and generate a 
nearly scale invariant primordial density perturbation, which has been tested by the observations 
of the cosmic microwave background radiation spectrum pj. The simpler examples are large 
field models of inflation with a renormalizable potential, a mass term plus a quartic interaction [2[ . 
When the quartic self-interaction dominates, the WMAP normalization of the primordial spectrum 
demands this coupling to be tiny, A ~ 10~ 14 . Present observational data also indicate that we 
live in an accelerated expanding Universe today Q, and that around 70% of the total energy 
density is made of a component, called dark energy, with negative equation of state w, close to 
-1 [H, El. This could be explained by a light rolling scalar field, usually called quintessence [1, Q. 
The dynamics of the field is such that at early times, during matter or radiation domination, its 
energy density pue is subdominant, and it is only today, when the field finds itself evolving in a 
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roughly constant potential, that pde dominates. Commonly, quintessence potentials are given by 
either exponentials or inverse powers of the field, i.e, by a non-normalizablc potential, which yields 
naturally a small classical mass of the quintessence field, below the Hubble parameter, without the 
need of fine-tuning. 

In trying to explain the origin of the inflation and quintessence potential from a particle 
physics model, it is unlikely that it will appear as an isolated entity with no interactions to any 
other degree of freedom. For example, inflation should be follow by a reheating period, such that 
we recover a radiation dominated universe at the end of it. The standard picture is that the inflaton 
couples to other light degrees of freedom into which it decays during reheating [7]. Moreover in 
warm inflation dynamics [8[ interaction of the inflaton field with other fields is needed to produce 
radiation concurrently with inflationary expansion. For models where the inflaton has sizable 
coupling to other fields, in order to maintain the required flatness of the potential, supersymmetry 
is typically used to control large quantum corrections Q. For quintessence models, there are viable 
models coupled to the dark matter fluid [Io| • Or quintessence fields can directly be coupled through 
standard renormalizable interactions to other bosons and fermions, providing for example time- 
varying masses for these degrees of freedom [ll|. But once coupled to other species, one should 
check the stability of the classical results against quantum corrections. 

These sources of quantum corrections have been already studied in the literature. For inflation 
for fields interacting sizably with the inflaton, their corrections usually must be controlled by 
supersymmetry. Often these effects are not a problem, but indeed they can either help with the 
inflationary trajectory, as in supersymmetric hybrid model [l2j], or lead to different observational 
signatures [l3|. For quintessence, in Ref. [l4| the one-loop corrections due to the self- interactions, 
among other cases, were studied by regularizing the theory with a cut-off A, and they conclude 
that for reasonable values of the high-energy cut-off, A < mp, quantum corrections do not spoil the 
quintessence potential. At early times, when the quintessence energy density is still subdominant, 
quantum corrections can be quite large compared to the tree-level potential, but still subdominant 
compared to the other components of the energy density. By the time the quintessence epoch starts, 
quantum corrections have already become negligible. A similar conclusion is reached in Ref. [15| 
for the self-couplings. More severe are the constraints on renormalizable couplings to fermions and 
scalars, which have to be really tiny and negligible to avoid distortion of the quintessence potential 

At first glance it seems that the previous studies on quantum corrections to the quintessence 
potential practically forbid couplings to other scalars and/or fermions. For inflation it appears 
that in order to have sizable interaction with other fields requires supersymmetry. Here we want to 
argue against these conclusions, and show that quantum corrections (at one and two loop orders) 
can be kept under control once the decoupling theorem is taken into account [HI . The result does 
not depend on whether we study a quintessence or inflation model. As such we will deal with a 
generic scalar potential V(<f>), where the field has a non vanishing background field value. We will 
focus on the interactions to scalars, by introducing another scalar \ with potential: 

^) = ^ + ^ 4 + ^V. (i) 

The potentially dangerous term is g 2 4> 2 x 2 > which can give rise to a large quantum contribution to 
the effective potential. This coupling induces a large field dependent mass for the x field, m x ~ g6, 
and when <p ~ mp , which is often the magnitude during inflation 0, [^, [l9l [2pj or quintessence [BJ , 
this can be larger than p 1 ^ 4 , p being the total energy density. However, if we do not have enough 
energy to excite such heavy states, physically we can expect them to decouple from the spectrum 
[l8l |. and their contribution to the effective potential to be highly suppressed. To properly address 
this issue, one must compute the effective potential using renormalised perturbation theory in 
which decoupling is already implemented. This is a well studied problem |2ll423l |. and we want to 
examine the implications of their results in the context of scalar field inflation and quintessence. 

The immediate problem in computing radiative corrections to the effective potential is that 
the latter must be renormalization scale independent; the point is how to choose that scale. One 
common approach is such that all log corrections due to different mass scales are kept small and can 
be resumed. For models with a single mass scale scale, the perturbative quantum corrections lead 
to logarithmic terms of the form ln n (m 2 / p 2 ) and the standard procedure for controlling the large 
log-terms is to choose the renormalization scale near the mass scale p *~ m. For a multi-mass case, 
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such as in Eq. ((T|), logarithmic corrections terms will arise of the form ~ In™ (m^/^t 2 ) In™ (m 2 //i 2 ) 
and if and m x are at very different scales, say m x ^> m^, there is no ideal choice of /i to control 
the large logs. However such terms arise in a mass- independent renormalization scheme, which is 
problematic for multi-mass cases for which there are disparate mass scales, since accounting for 
decoupling effects can not be done. It is physically better motivated to use a mass dependent 
renormalization scheme [HI , in which any field with a mass much bigger than others in the system 
has its quantum effects suppressed in powers of the light-to-heavy mass ratio. Thus in a mass 
dependent scheme, quantum corrections lead to terms in the perturbative expansion moderated by 
power suppression, ~ (ji/m x ) ln n (m 2 i //i 2 ) ln m (m 2 /^ 2 ). In this case the choice of renormalization 
scale /i 2 -C rn x , although will lead to large log terms, it is not important since the power suppression 
term dominates the perturbation series and keeps it under control. 

Resummation of the logs is done by app lying renormalization group (RG) techniques to obtain 
the RG- improved effective potential [25|-|30|. For multimass scale problems, the prescription given 
in Refs. [2l|, |23j is to always choose the renormalization scale to be of order the lowest mass scale, 
because the decoupling theorem will ensure that the heavy mass states do not contribute. This 
can be seen directly when using a mass dependent renormalization (MDR) scheme, because at any 
order the logs are modulated by the appropriate threshold functions [22|, |24| . Recall that these 
threshold functions suppress those massive contributions that are much above the renormalization 
scale. Therefore, in the MDR scheme one can immediately see that one gets the most rapid 
convergence of the perturbative expansion by choosing the RG scale below all thresholds. 

Still, the parameters in the effective potential have to be specified at some renormalization 
scale, and there is no direct information in the potential about how to choose that scale. A 
large value of the potential, or its curvature, etc., has no direct relation with the choice of the 
renormalization scale. One needs outside information on the model. For example, when computing 
the Higgs effective potential one does not know the self coupling parameter at the electroweak scale. 
One procedure then is to invoke a higher symmetry, like a Grand Unification Theory (GUT) model, 
which then specifies the coupling at the GUT scale, and use the renormalization group equations 
(RGE) to run it down. This is the standard procedure when studying the effective potential in both 
the Standard Model and the Minimal Supersymmetric Standard Model. If there is no theoretical 
argument from the symmetries of the model, another possibility would be to do some scattering 
experiment from which to extract the value of the self-coupling at some momentum scale, similarly 
to what is done for gauge couplings 3l| • Then, repeating the experiment at different energy scales 
one can further confirm the predicted running 31 1. 

In inflation and quintessence models most often they are not embedded in a higher theory, 
and moreover outside phenomenological information about the parameters is unavailable. For 
example, in chaotic inflation normalising the potential to the cosmic microwave background (CMB) 
amplitude leads to a tiny self-coupling of A^ ~ 10 -14 [l|,|2|]. This was all done through a classical 
calculation. What renormalization scale has this specification of A^ ~ 10~ 14 been done at? No 
information internal to this calculation tell us about the scale. This is a common problem in this 
sort of inflationary model building. As we will discuss in this paper, this implies an arbitrariness 
in the model predictions. 

In the case of our model, there are three different couplings that need to be specified, the <fi 
self-coupling A^, the x self-coupling A x , and the coupling g between the two scalar fields. If we 
could scatter particles, that would give the couplings and masses at that scale, but that possibility 
is not available. Indeed, this problem is general to any scalar field inflation and/or quintessence 
calculation. One has no a priori information about the values of the couplings at any given 
renormalization scale, so it is a matter of choice, which leads to considerably ambiguity in the 
predictions from any given model. To highlight this point about the ambiguity, decoupling has 
an interesting implication. If one chooses all the parameters at a scale well below all thresholds 
in the model, then there will be no quantum correction to the effective potential as such, and 
thus the tree-level potential is the full effective potential. This is the implication of decoupling in 
application to the effective potential in Refs. |21H23| . On the other hand, if there are massless 
states, or one only knows the value of the parameters at some scale intermediate to the masses in 
the theory, then one must implement the renormalization group improved potential. 

In our model, we have three couplings A^, A x , and g, and the only thing we know is that 
either for inflation or quintessence must be very tiny, although A x and g are unconstrained. We 
also know that the x mass is heavy when the inflaton background field value is large, whereas the 
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inflaton mass must be smaller than the Hubble rate, or in the case of quintessence field, it must 
be very tiny. Following the above logic of a scattering experiment, the energy scale certainly will 
be below the Planck mass, and an internal consistency for our model would require to be very 
tiny. Then decoupling requires that heavy states yield no radiative corrections. 

In this paper the decoupling concept will be developed and applied to inflation and 
quintessence models. In section [TT] we will first discuss quantum corrections due to renormaliz- 
able interactions to other scalar field, computing the RG improved effective potential at one-loop 
in a MDR scheme. Details about the calculation of the RGEs in the MDR scheme are provided 
in Appendix A with one-loop results, Appendix B with some two-loop results and the one-loop 
effective potential is given in Appendix C. In section HTT1 we discuss the non-renormalizable self- 
interactions relevant for a quintessence field. In order to check the stability of the potential against 
quantum corrections, we study in section IIVI an example for an inflation model example, and a 
quintessence potential. We present the summary in section Ivl 



II. ONE-LOOP CORRECTIONS: RENORMALIZABLE INTERACTIONS 



In order to proceed, we will first check the one-loop corrections due to the renormalizable 
interactions of an extra x field coupled to the scalar field <f>. We assume that <j> is a light field with 
a non- vanishing vev, whereas the x state is heavy (heavier than the Hubble expansion rate), with 
zero vev. The tree level potential is given by: 

F(°) (0, x) = + \m\<? + ^ + \m 2 x X 2 + ^X* + ^ 2 xV , (2) 

Loop corrections can give rise to a cosmological constant term f2, a quartic self-interaction A^ and 
mass term m^, and for consistency we include them already at tree- level. The aim is to show that 
those interactions do not pick up large corrections due to the heavy field, and they can remain as 
small as required during the inflationary or quintessence phase. We do not couple the <f> directly 
to fermions for simplicity, and focus only on scalar couplings, but we allow a Yukawa coupling of 
the x field to Np massless Dirac fermions, 

£vuk = -hxipip- (3) 
One- loop corrections are given by the Coleman- Weinberg potential 32] : 

A^ (1) = i fAYMf + rt), (4) 



2 J (2tt) 4 ^ 

where a — 0, x, and M a are the field-dependent masses: 

M 2 X = g 2 4> 2 + m 2 x , (5) 
Ml = ^ + ml , (6) 

(We will denote <j> for both the quantum field as in Eq. (2) and in all the Lagrangians in this paper 
as well as for the background field value as in the above expressions for M x and M^, since the 
correct usage will be obvious in each case.) The divergent integrals in Eq. (j4]) can be regularized 
by using a cut-off A, and keeping only terms that do not vanish when the cut-off goes to infinity 
we have: 

AVU=Y(^ 2 + ^(^-~)) (7) 

a 

For renormalizable tree-level potentials, the cut-off divergent terms are subtracted by adding coun- 
terterms: 

V c M)=^ + \5ml^+ 5 ^ (8) 
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5 ^ = ^Z2( 3A ^ + 12 5 4 ) ( ln — ~ 1) ) ■ (11) 



and imposing suitable renormalization conditions at some arbitrary scale y, on the effective potential 
[32| . This allows to remove the quadratic and logarithmic divergent term, by choosing: 

5Q = -^K + -x) + ^K + -x)( ln f -i)) . (9) 

5ml = + ^ + 3^2 ( A * m * + 2 9 2m D ( ln ^ - 1)) ' ( 10 ) 

3^ + 12/) (ln| 

Adding T4t to Eq. (|7|), the divergent terms cancel out and one is left with the logarithmic contri- 
butions depending now on the renormalization scale \x\ 

^ m = ^5>«('»f -£)• a*) 

Notice however that there is an arbitrariness in choosing the finite terms in the counterterms, and 
different renormalization conditions lead to different finite contributions [331 ] also in the effective 
potential. The above prescription has been chosen in order to match the standard result for 
the 1-loop effective potential using dimensional regularization and minimal subtraction (MS) as a 
renormalization prescription. 

Independently of the implemented renormalization scheme, physics cannot depend on the 
arbitrary renormalization scale //, and the effective potential V = + AV^ has to satisfy the 
renormalization group equation (RGE): 

w =(4 + ^^-^4- 7xx l) y=o ' (13) 

where A a denotes both renormalizable couplings and mass parameters in the potential, j3 a = 
3A a /<91n/i their beta functions, and 7^, 7 X are the anomalous dimensions of the fields 1 . 

The solution to the RGE then provides the RG-improved effective potential [2l|, |25U30| , given 

by: 

V{<t>, x, A Q ; m) = V(4>(t), X (t), A«(t); e*n) . (14) 

We can now evaluate the effective potential at any given scale t by appropriately changing fields, 
couplings and masses: <fi(t), x(t)i -HO are now running parameters, with scale dependence (t- 
dependence) given by the corresponding RGEs. As remarked in [26], order by order in perturbation 
theory, the Lt/i-to-leading log order RGE improved effective potential is given by the (L + l)doop 
RGE functions, and the L-loop effective potential at some boundary value of t. The main idea of 
the RG improved method is that by choosing adequately t, the potentially large logs appearing on 
the LHS of Eq. (fT4)) can be resummed. Thus, at lowest order the RG improved effective potential 
reduces to the tree-level potential with couplings, masses and fields given by the 1-loop running 
parameters: 



= Q(i) + i m 3(t)^) + ^^( f ) 

(15) 



l~,2/j.\...2/V\ 1 ^tW„4/j\ 1 1 „2/jA, 2/j\ il 



m 2 x (t) x 2 (t) + -^x\t) + -g 2 (t) x 2 (t)4> 2 (t)^ _ 

The key point is how to choose t*, i.e, which is the best choice to evaluate the effective potential. 
As stressed in Ref. [23[ , one can choose either a value at which the one- loop potential has the least 
/i-dependence, i.e. 

D(V {0) + AV {1) ) =0, (16) 



The anomalous dimensions are the logarithmic /i derivatives of the wave-function renormalization constants of the 
fields . We have included in Eq. I|13| l the contribution of -y x for the sake of generality, but we only consider 
situations when \ = 0, X referring to the vev of the field. 
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or the scale at which the loop expansion has the best apparent behavior, i.e., AVW(i*) = 0. The 
optimal situation occurs when both criteria are met by the same choice t*. This can be done in a 
model with a single mass scale, say a scalar field <j> with self-interactions and no couplings to other 



fields, so that the choice (25l, [2Q. I28l430l : 



U = hx^. = ±U-^-, (17) 



fulfills both conditions, and is equivalent to evaluating the one-loop potential at the scale ijl = M^. 
But in the presence of very different mass scales, say <C M x , the choice is not obvious. The 
problem is how to rearrange the loop expansion in terms of small parameters for the series expansion 
to make sense (and to be resummed). With two different mass scales, the correction to the effective 
potential can be written as [2l|, l22| : 

A ^ = .^ 4 E iLi )T. F ^MiK^/9 2 ]s^, (is) 



l=L-{i+j) v ' i,j 



where 



q 2 M 2 



By taking as a boundary condition either — or s x — 0, to evaluate the L-loop effective 
potential at given order, there are still potentially large contributions due to \nM 2 /M^ in Eq. 
(1181) . and the series cannot be truncated at any order. 

However, as remarked before, one does not (and should not) expect heavy states to modify the 
low energy physics, and thus the main issue to address is how to incorporate decoupling of heavy 
states in the improved effective potential [21-23], that is, how to get rid of the troublesome large 
logs by adopting a physical condition. For example, in Eq. (|15p the decoupling can be incorporated 
in the running parameters through their RGEs. In a mass-independent renormalization scheme, 
with no reference to mass scales, decoupling has to be implement by hand by the use of step 
functions in the RGEs and matching conditions for masses and couplings at each threshold: for a 
given state, when its mass becomes heavier than the renormalization scale, its contribution drops 
from the RGEs and the effective potential. In particular we will have [23] ]: 

AV« = ]T^ 2 -M2)K«, (20) 

oc 

64tt 2 \ n 2 2 J ' 1 ' 

and the optimal choice for i* (or /x*) is then given by the lower threshold in the model. Indeed 
in that case all massive states are decoupled, so that Af' 1 '^,) = 0, i.e., the effective potential 
is given by the tree-level potential with masses and couplings evaluated at low energy. The effect 
of the heavy states appear when integrating the RGEs from high energy down to the low energy 
regime. 

However, threshold effects are more naturally taken into account when adopting instead a 
mass dependent renormalization (MDR) scheme (22J. Following Ref. [2~3 |. effective couplings and 
masses are defined after subtracting the divergences (regularized using dimensional regularization 
for example) of the 1PI Green functions, by imposing suitable normalization conditions at the 
euclidean external momentum p 2 E = —p 2 = /i 2 , with ji 2 being the arbitrary renormalization scale. 
The Appendices give a detailed account of the MDR scheme, beta functions and 1-loop correction 
to the effective potential in this scheme, while the following just summarizes the main results. 

We have already mentioned the scheme-dependence of the finite contributions in the countert- 
erms, or equivalently the renormalization constants Z a [33|: 

6m 2 a = mKz-l-1), (22) 
SX a = X a (Z Xa -l), (23) 

sn = n(z n 1 - 1) , (24) 
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that relate the bare parameters (denote by a subindex "0") with the renormalized ones: 

m a(^) = Z <t> Z rnim 2 aQ , (25) 

\ a (ji) = Z$Z£\a, (26) 
= Z Q Ct . (27) 

By modifying the subtraction conditions, we are explicitly including finite contributions from the 
1PI functions, which carry the dependence on the different mass scales of the model. When taking 
the derivative with respect to the arbitrary renormalization scale /i, that dependence appears 
in the beta functions as threshold functions depending on the different ratios M 2 /fj, 2 . These 
threshold functions modulate the contribution of each massive state to the running of the different 
parameters. And in the exact decoupling limit /i = the threshold functions vanish. For the 
renormalizable parameters in Eq. ([2]), the beta functions for the couplings and mass parameters 
are given by: 

(4tt)X = 3A^ 2 (a ) + 12 5 4 ^ 2 (a x ), (28) 
(4tt) 2 /3 Ax = 3\ 2 x F 2 (a x ) + l2g 4 F 2 (a 4> ) + (8A x /i 2 - 48/j 4 )7V f , (29) 
(4^) 2 /3 3 2 = g 2 (X 4> F 2 (a 4> ) + X x F 2 (a x )+8g 2 F 1 (a x )+4N F h 2 ), (30) 
(4n) 2 l3 h 2 = h 4 (6F 3 (a x ) + AN F ), (31) 
(4^) 2 /3 m 2 = 4N F h 2 m 2 x , (32) 
(4^) 2 /3 m 2 = 0, (33) 
(47r) 2 /fe = 0, (34) 

where Np is the number of massless fermions, and we have defined a a — M 2 /[i 2 . Notice that 
by using a MDR scheme, with no couplings of the field <j> to fermions, the mass parameter 
and the vacuum energy contribution il do not run: they are fixed by the boundary conditions. 
That is, the pure (quadratically) divergent terms from vacuum diagrams with no reference to the 
external energy scale can be subtracted from the bare parameters in the potential, leaving a fixed, 
finite contribution. For example we can always impose as normalization condition for the vacuum 
contribution: 

f2(^)=0. (35) 

On the other hand the effective field dependent mass M 2 does run, due to the running of the 
coupling constant A^. 

The expressions for the threshold functions Fi, i — 1,2,3, are given in Eqs. (|99| - (1101[) 
in Appendix A (also some two- loop expressions are given in Appendix B). They can be well 
approximated by: 

FM^ ,' + C T 2 , (36) 

with the coefficients c^- for each function given in Table HU In the massless limit, all threshold 
functions reduce to one, which recovers for the couplings (dimensionless parameters) the standard 
RGEs computed for example in MS (see Appendix A). On the other hand, when the ratio a a ^> 1, 
we have Fj(a a ) ~ 0(l/a a ), i.e., power suppression of the heavy state contribution. Decoupling 
is not instantaneous, as can be seen in Fig. [T] Threshold functions smoothly interpolate between 
the high energy regime where massive states can be viewed as massless, and the low energy theory 
without heavy states. 
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FIG. 1: Threshold functions: F 1 (m 2 /p 2 ) (solid, black), F 2 (m 2 /p 2 ) (dashed, red), and F s {m 2 /p 2 ) (dot- 
dashed, green). 



All that remains now is to fix the initial conditions at some scale fj, to integrate the RGEs and 
obtain the values of parameters when all masses decouple. With those we evaluate the tree-level 
potential to obtain the RG-improved effective potential at 1-loop. It is shown in Appendix C that 
by substituting back the solution for the running couplings at low energy in the tree-level potential 
one recovers the one-loop correction computed in the MDR scheme: 

i / \ 2 (h A M 2 M 2 \ 

= -fi ~ HM 2 ,, 2 )) + 5 V(ln -± - I(M^))J , (37) 

where 

1(a) = In a - 2 - yT+laln . . (38) 

Vl + 4a + 1 

The main difference between the one-loop correction computed in a mass independent renormal- 
ization procedure, and the MDR scheme one, comes into the non-logarithmic contribution, and 
it is due to the different scheme-dependent finite contributions in the renormalization conditions. 
Comparing Eq. (|3"?| with Eq. (fT2")l . the constant "3/2" term is replaced by a threshold function 
1(a), which controls the contribution of the original log term. Thus, whatever the hierarchy among 
the masses, we obtain the exact result at 1-loop: 

AV (1) (n = 0) = 0. (39) 

In section IIVI we will present some examples of the procedure for an inflation model, and 
a quintessence one. We want to check the impact of the radiative corrections on the infla- 
ton/quintessence potential as 4> changes. Notice that by changing cf> we are implicitly changing the 
threshold conditions that depend on the effective masses, and therefore the values of the couplings, 
so effectively what we have are field dependent couplings. We will check that at least during the 
regimes of interest, when </> evolves the couplings remain in the perturbative regime and none of 
them picks up a large correction. For example for the quintessence model, we can impose that we 
have indeed a quintessence regime, such that for <fi 3> rap the quartic coupling is tiny (or even 
cero). By going backwards in time, i.e., taking smaller values of the field, the evolution should be 
such that the coupling never gets large enough to disturb the standard evolution of the quintessence 
field. 
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III. ONE-LOOP CORRECTIONS: NON-RENORMALIZABLE SELF-INTERACTIONS 



Most Quintessence potentials (and some inflation potentials) are typically given by non- 
renormalizable (f> potentials, Vath (</>), so that the tree- level potential is now given by: 

V (0 H^x) = V NR (t) + \m\<? + ^ + \m\^ + ^ + \g\^ , (40) 

The ^-dependent mass of the <p field picks up an extra term due to the non-renormalizable inter- 
action, and in Eq. (O we have to replace Mj by V^ R + Mj. The 1-loop correction AV^ 1 ) can 
be split into a renormalizable and a non-renormalizable contribution, owing to the origin of the 
field-dependent masses: 

AW 1 )=AyW+AyW, (41) 
where AV T e2 is given in Eq. ([?}, with M a = M^, M x , while AV^p is given by: 

Vnr ~ 32tt 2 A + 64^ [ V nr^ nr + 2M 4) ) ^ln A2 ^+^ln ^ j(42) 

Having dealt with the renormalizable interactions in the previous section, we come back to the 
quadratic cut-off and log dependent term due to the non-renormalizable interaction. The standard 
approach, which we will follow, is then to consider A as the effective ultraviolet cutoff for the 
model, such that: 



r(i) ^ Vnr a 2 



AV W _ A 2 (43 s ) 



where in AV^^ we have kept only the dominant quadratic contribution |16j . 

Without loss of generality, let us consider a generic quintessence potential of the form: 

The scale M is model dependent, and we do not consider any particular value; here it only 
parametrizes the value of the field well outside the quintessence phase, <j) ~ while the 
quintessence regime happens for values of the field <f> 3> mp. Now we want to check that in- 
deed AV^, does not provide large corrections to the potential, AV^ <C Vnr- This condition is 
not difficult to fulfill when <fr ^> M, and in this regime we simple have AVmr/Vnr — (A/</>) 2 , so in 
the quintessence regime when <f> > mp, the mass squared V^ R has become tiny, and the effective 
ultraviolet cut-off can be taken close to the Planck scale and still A/0 < 1. However, at early times 
when <f) ~ M <C mp we have that AV£l/V NR ~ (A/M) 2 , which can be large unless the cut-off is 
well below M. 

Nonetheless, one can argue as done in Ref. [l4[ that this is all right as far as the one-loop 
contribution is suppressed not with respect to the tree-level potential, but with respect to the 
dominant energy density at the time. For example there are some restrictions on the amount of 
dark energy at the time of BBN, with pde < 0.2p ra d, which in this case should be satisfied by 
the 1-loop effective potential. At earlier times, the quintessence field may find itself fast rolling 
the potential, with the Universe dominated by its kinetic energy density pke (kination). The 
condition to be on the safe side would be that still AV^ R <C Pke- However, it is not clear how 
to reconcile a fast rolling field with the calculation of the improved effective potential, and the 
approximation may break down. Because of that we do not pursue the calculation of the effective 
potential into that regime. Whenever a kination phase due to the quintessence field in the early 
universe, we can check that after that quantum corrections do not mess up the evolution of the 
quintessence field, and that quintessence domination is reached today. Going backwards in time, 
if the quantum corrections are subdominant by the time of kination, we assume that they will not 
grow as much as to change this phase. We have no means to consistently check this assumption, 
but we consider it a reasonable working hypothesis. 
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IV. RESULTS FOR INFLATION AND QUINTESSENCE 

A. Inflation 

Let us consider inflationary potentials of the form: 

V(d>) = X^(-^J , (45) 

which is added to V(x) in Eq. (JXJ) . Due to the coupling between <p and x, radiative corrections 
will always induce a quartic interaction for the inflaton field. Just to keep the discussion simple, 
we focus on the n = 4 case. The tree- level potential reduces then to Eq. @, with 2 = 

In such chaotic potentials inflation takes place for (j) > mp, and therefore the x held gets a 
large mass g<p — mp, for moderate values of the coupling g. On the other hand, the inflaton mass 
is M? = < H 2 <C M 2 . . Then, following our prescription, the appropriate renormalization scale 
fj, = H* for examining physics during inflation is of the order of M$. At this scale, the threshold 
functions imply a suppression of the effect of the \ loops in the renormalization group equations. 

Notice that in the standard MS scheme, the generic approach is to suppress large logs, in which 
case at one- loop order it would mean taking (j,* ~ M x . However, in the MDR scheme the large 
log contribution is suppressed by the threshold function prefactor when fi is chosen at the lowest 
threshold. In fact, at two- loops at order g 4 , there will be large logs depending on both masses 
and M x , and in the MS scheme it is not clear which is the best choice for /z*. Nevertheless, the 
MDR scheme is unambiguous that it is around the lowest mass scale. 

Up to order of magnitude, the above approach fixes the choice of fi* in the MDR scheme. But 
there is an uncertainty in the exact value one should choose. This underlies an inherent ambiguity 
in the value of the effective potential, which ultimately implies a theoretical uncertainty in the 
coupling and thus on the model predictions such as the amplitude of primordial perturbations. 
If the effective potential could be computed exactly, then it would be completely /i-independent 
and any choice of /i* would be equally good. However in any real calculation, where the effective 
potential is calculated only to some finite order, often just one-loop order, the choice of /j,* must be 
made carefully. For a given choice of /j,*, slightly larger or smaller values should result in the same 
answer, and if they do not, the selection of /i* is flawed for the given order in the loop expansion. 

To implement the calculation of the RG improved effective potential, two renormalization 
scales in general are needed. First is the scale, which will be called p, where the initial values 
of the parameters in the theory are specified and second is the scale where we want to use the 
potential to do physics, which we have already denoted /i*. A further detail in specifying the initial 
values of the parameters at p is this specification in general must be given over some range of (f>, 
and thus the masses in the system. As <j> changes, these ^-dependent masses will change, thus this 
choice of p and/or initial values of the system parameters can change. 

There are two approaches we will consider for initializing the RG improved calculation of the 
effective potential, which will be referred to as the high and low energy approaches. In the low- 
energy approach the parameters of the system are initialized at the scale where one is interested 
in using the effective potential to do physics, thus p — fi*. Since this /i* in MDR scheme is at the 
scale of the smallest mass, it means there would be very small quantum corrections to the effective 
potential from heavy mass states. Moreover as <p changes, thus the masses in the system, in 
principle one could use RG to move the value of /z* to optimize the quantum corrections, although 
the effects would be small. In the high energy approach, the values of the initial parameters are 
specified at some high renormalization scale p and then the renormalization group is used to run 
the parameters down to the scale n*, where physics is to be done. 

An example of implementing the low-energy approach in the case of inflation would be to fix 
the value of at the epoch of inflation corresponding today to the largest observable scale. This 



2 We have shown in section [II] that the mass parameter does not run in the MDR scheme, so that imposing as 
a boundary condition m^,{fi) = ensures that this parameter vanishes at 1-loop at any other energy scale. 
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FIG. 2: The running of A^, with g 2 = 10~ 4 . The red star indicates p = p; at this scale some assumed 
physics informs us that A^> is constant over the considered range in (f>. Blue dots indicate the position on 
each curve where fi — /i*, and where the parameter should be taken to improve the effective potential. 
Curves with smaller <f> than those shown quickly drive A^ negative before n* is reached. At larger <j> the 
curves are flat and constant, with a value that increases only slightly as <f> increases. 

value of A^ can be determined from density perturbation constraints from measurement of the 
cosmic microwave background. The renormalisation scale this corresponds to would be (following 
our prescription) the lowest mass scale in the theory, fi* ~ M^. We may therefore say that at 
these values of 4> and /i, observations tell us A^(/i*) s» 10~ 14 , and this is all. If we have a more 
complicated model with additional parameters, we must also be able to specify the remaining 
parameters such as g 2 also at /i*. In this approach there would be no or very small quantum 
corrections, so that the tree-level potential would be almost identical to the RG-improved one. 

The high energy approach might be implemented if the scalar potential were embedded in a 
higher theory, and some symmetry at a high energy scale p specified the value of the parameters 
and over some range of <j). Then the RGE could be used to run the parameters down from p to /i* 
where one wishes to do physics with the effective potential. 

An important point here is the matter of initial conditions is not simply a mathematical 
concern. There is in general missing physics in inflationary models. The predictions one obtains 
from such models depend on this missing physics. Thus for a given inflationary potential, depending 
what higher theory it is embedded in, different specifications might emerge for the value of the 
coupling A^ at some high energy scale, which when evolved down to /i* will lead to different 
predictions for inflation and large scale structure. This can be viewed in an alternative fashion. In 
the low energy approach, we require the inflation model to be consistent with observation, thus the 
parameters are determined by the physics of inflation. This might then be used to place constraints 
on the unspecified physics at higher energies. 

Let us now examine the behavior of the parameters with RG running. For the high energy 
approach at some renormalisation scale fi = p some unknown physics specifies that the parameters 
(and A^ in particular) maintain their values over a range in 4>. Although we are not considering 
any specific model for the high scale physics, we would like to investigate the procedure for how 
this would in principle be done. As such, what we will do to determine the values of the parameters 
at scale p is run from /i* (where the parameters are specified or known) to p, at the value of <fi 
where a constraint on the parameters exists. For our purposes this will be <j> = m p . This is just 
to ensure that the effective potential will match observed constraints at <f> = m p . We would like to 
check it continues to do so at larger/smaller values of <fi. 

Thus we start with the following parameters at = m p and fi = (j,* = M^: = 10 -14 , 
g 2 = 1(T 4 , A x = 1CT 3 , h 2 = 1CT 4 , m\ = l(r 9 mf,, m$ = 0, Np — 8, and p = lCT 5 m P . We run 
the parameters upwards to p and find the value of the parameters there. At this scale our assumed 
physics keeps these parameters unchanged with <j). Thus we now change 0, and run back down to 
p = [i*. This lets us probe the effective potential at different values of <fr. 

What is found is most parameters remain nearly constant (varying less than 0.1% over the 
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FIG. 3: Curves showing how A^/i*) varies with <f>, for different values of fi. As before, g 2 = 10 -4 . Since 
A^,(/x*) appears in the RG- improved effective potential, this is the physically important quantity. Curves 
with fl larger than those shown quickly drive A^(/i*) negative before <f> is as small as 0.1m p , and have 
mildly larger values at large <j>. At larger (f> the curves asymptote to a constant value. Increasing g 2 has a 
similar effect to increasing fi 2 , and the increase is one may be compensated for by a decrease in the other, 
allowing for larger values of g 2 than plotted here. 

range in /i considered); the only parameter that runs appreciably is A^, and is plotted on Fig. [5] 
The blue dashed line on Fig. [2] shows this initial curve (</> = m p ) for A^. The red star denotes p., 
which is at the same location on all curves. 

Recall at this scale fi, the parameters are assumed to be independent of cj>: through this 
assumption we are including the missing physics. We run down from the red star to the low 
scale fi* for different values of <f> (the blue dots denote the point ft,* on each line). At this scale, 
the parameter should be taken and inserted into the tree-level potential to generate the one-loop 
improved effective potential. 

It is clear the curves are very flat at the scale /it*; it makes little difference if the parameters 
are evaluated at precisely fi* or within some order of magnitude of this scale (or indeed, many 
orders of magnitude below it). This is precisely as we would expect, as we have chosen fi* just for 
this property, so as to satisfy the /i-independence of our improved effective potential. 

As each curve generates a single value of \(fi*) for a single value of 0, we can construct curves 
of how A(/i*) varies over <f> directly. To do so we compare the values at ft* (the blue dots) for a 
large number of curves. Furthermore, we can examine how these curves change if the value of ft 
(the horizontal position of the red star) is changed. The result is Fig. [31 maintaining the above 
parameters but varying fi. 

How can we understand this behaviour? Consider the equation governing the running of A^, 

Eq. d2HD, 



(47r) 2 /3 A , = 3A^ 2 (a ) + 12 3 4 F 2 (a, 



(46) 



We will now make some approximations. If we assume fi 2 <C M 2 (certainly true for the numerical 
range shown above), we may use the approximate behaviour of the threshold function, F 2 {a x ) m 
fi 2 1 \g 2 (f> 2 ) . Despite this suppression, the value of A^ is still small enough for the second driving 
term to dominate and we find: 



,dA 



(47) 



We can solve this if we make the assumption g 2 is constant, which is well supported numerically. 
The solution is 



Mm) = Mm) 



6g 2 



(4tt) 2 \(j> 



(48) 
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We see now precisely the behaviour in Fig. ([2]); starting from an initial value at A^,(/2), A^ quickly 
reaches a constant value when fi -C p. We also see the </> dependence exhibited in Fig. As 
(j> is steadily reduced, the second term becomes larger, eventually dominating the first term and 
driving A^ negative. Larger values of p have the same effect as this. For A^(/i*) to be stable 
against changes in (f>, we require the second term smaller than the first. This gives the constraint: 



6g 2 ( ^ 1 



(47T) 2 



« A (A) , (49) 



which can be written as: 



3/2 > 10M 4 {p) . (50) 
For inflation, with \<p(p) ~ 10~ 14 and at its smallest <f> — 0.1m p , we find: 

3/2<10 12 GeV. (51) 

This agrees well with Fig. [3] if g 2 s» 10~ 4 , then p may be no larger than 10 14 GeV. 

When inflation ends at around 4> ~ m p , the suppression is reduced compared to earlier, larger 
field values. As the field strength continues to decrease, the effects become gradually larger. The 
values of p and g 2 in the plots have been chosen to display this behaviour between <j> = m p and 
<f> = O.lrrip. 

For completeness, note that if a large mass parameter is given to x, such that M x m x ^> 
g 2 <j> 2 , the solution becomes independent of (j> at these field values. However, for the solution to be 
valid, we would require M x m m x p to suppress the x loops. 

Let us now pause and consider the physical interpretation of these results. Without decoupling, 
we would expect corrections to A^ of size 0(g A ). Instead, due to the suppression of the x-loops, we 
find this result is reduced by a factor of p 2 /M 2 . Deriving an effective field theory, by integrating 
out the heavy degrees of freedom from a more fundamental theory, typically leaves behind a low 
energy theory along with terms in the potential suppressed by factors of O (E 2 /m 2 ), where E is 
the low energy scale the theory is probed at, while m is the energy scale associated with the heavy 
degrees of freedom. For instance, classical electromagnetism encounters corrections with a size 
O (E 2 /m 2 ) where m e is the mass of the electron. We might have expected in the case of inflation 
to find quantum corrections suppressed in a similar fashion, with corrections O (E 2 /M^j. But 
when dealing with an effective potential directly, it is not clear what energy scale to associate to 
E. The effective potential is computed from diagrams with external legs set to zero momentum, 
and the only other scale is the renormalisation scale, upon which our results cannot depend. We 
now have an answer (at least for the form we have assumed our missing physics take): it is the 
scale associated with p, where the theory is set. If the extremely large field-strength associated 
with the inflaton (with <f> > m p in chaotic inflation) serve to generate field dependent mass of the 
X-field far above this scale, its contribution to the effective potential is suppressed. 

At what scale should we expect p to be? Inflation is associated with high vacuum energy 
densities. But the relationship between the renormalisation scale and physical energy scale is not 
clear. Without an obvious way of relating these two properties, it is difficult to motivate any 
particular choice of p above any other. 



B. Quintessence 



Now we revert to the generic quintessence potential introduced in Eq. ([4"4"|). Quintessence 
occurs when </> 3> mp, and thus, any field at least moderately coupled to <ft wn l acquire a large 
mass from the <fi background field value. Therefore, again one expects only power corrections to 
the effective potential from such heavy states. From the previous example, we can perhaps already 
guess that we should expect the typical g 4 corrections to A^ to be suppressed by 0(p 2 /M 2 ). 

In a similar fashion to inflation, we have only limited knowledge of the effective potential. 
We know only that today, with 4> ~ m p , the value of A^ appearing in the effective potential 
must be small enough so that the dominant contribution comes from the non-renormalizable term 
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FIG. 4: The running of A^, with g 2 — 10 -20 . As with inflation, the red star indicates p where the 
parameters are assumed to be <f> independent. The blue dots on each curve indicate /x* for each value of 
(p. Increasing g 2 has a similar effect to decreasing <f> 2 . 



VnrO 1 )- Another way of stating this is that the tiny effective field mass of ~ 10 33 eV remains 
unchanged by the size of (//*). From this, we can place a constraint on X^(fi*) <C 10 -124 . 

At earlier epochs in the history of the universe, the quintessence field strength is smaller. As 
discussed earlier, at very small field strengths the effective potential description breaks down, but 
we would like to make sure the potential is not disrupted by quantum corrections over at least an 
order of magnitude in (f>. 

We proceed in just the same way as in the previous section with inflation. The system of 
equations are initialised at \i = /i* with <f> = m p . We take parameters at this scale as: A^ = 10~ 126 , 
g 2 = lCT 20 , A x = 1CT 3 , h 2 = 1(T 4 , m 2 = l(r 9 m 2 , m 2 = 0, Np = 8, and p = lCr 57 m p . All is 

identical as with inflation, except for the size of A^, g 2 and p, and of course the non-renormalizable 
contribution to M<f,, The result is Fig. [H where once again the blue dashed line indicates the 
initial curve (4> = m p ) that we use to find appropriate parameters at p, so that (by construction) 
we match observational constraints at this value of cj). The red star indicates p where the initial 
curve ends, and all other curves (each for a different value of <f>) begin. Blue dots indicate the 
location of /x* on each curve. 

The results are very similar to inflation. Smaller values of 4> cause A^/z*) to rapidly approach 
zero. Due to being proportional to an inverse power of <p, decreasing 4> also increases /i* = M^. 
Thus the blue dots on each curve shift to the right as (j> is decreased, the opposite for inflation. 

Making 4> small enough eventually drives A^,(/i*) past zero. Our concern for quintessence 
is different to that of inflation; we are not concerned with A^ becoming negative; rather we are 
concerned with it becoming large compared with the non-renormalizable term in the potential. 
However, once A^/it*) becomes negative, smaller values of 4> decrease it further, such that its 
absolute value becomes large enough to be problematic. Thus, to ensure the viability of the 
quintessence model, it is sufficient to make the same demand as with inflation: that A^, remains 
positive for a a considered range in ip. 

Just as before, we construct curves of A^(/i*) with respect to 4>, shown on Fig. 0J As with 
inflation, very small values of <fi likely signal a breakdown in the description of the field in terms 
of an effective potential. This is particularly clear in the case of quintessence, where in the early 
universe and at small field values (f> m p the field may be moving extremely quickly, to the 
extent that the field's behaviour is dominated by the kinetic terms. Again at larger field values, 
there is little difference to the value of the parameter. Due to the non-renormalizable term in the 
potential, smaller values of <f> actually leads to a larger value of /j*. This makes little difference to 
the behaviour of A0(/i*) with respect to <p however, as the curves are so flat at low values of /i. 
Fig. [5] shows this dependence of A<^(/Lt*) with <j), for different values of p. The behaviour is much 
the same, and can be understood in precisely the same way. 

The equations remain unchanged when moving from inflation to quintessence. The solution 
is therefore identical, and following the same line of reasoning, we may therefore write down the 
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FIG. 5: Curves of X t j > (fj 1 *) against p. The value of g 2 remains at 10~ 20 . At large values of <p the curves 
become flat, while small values induce large changes in A^(/i*), eventually driving it negative (and even- 
tually, large). As with inflation, larger values of g 2 may be compensated for by using smaller values of p, 
and vice- versa. 



constraint: 

gp < WM^p) ~ 10~ 42 GeV . (52) 

Exactly how constrained the coupling g 2 is, depends upon the renormalisation scale p at which 
the physics is set. 

Quintessence is a low energy phenomena, operating at scales great many orders of magnitude 
lower than inflation. The constraint on means that the combination gp must be much smaller 
than in the case of inflation, to avoid generating corrections to the effective potential that ruins 
the tree-level behaviour. However once again there is a great deal of ambiguity as to what an 
appropriate value of p should be. This is information not included in quintessence models, just as 
it is not included in models of inflation. 



V. SUMMARY 



The approach in this paper is consistent with the concepts of low-energy effective field theories 
[35M37| . In the effective field theory approach, at low energies, reflected by a low renormalization 
scale, the effective theory is obtained from the full theory by removing all heavy mass field inter- 
nal propagator lines and treating the heavy fields only as external particles in Green's functions, 
with all other effects emerging through renomalization of the parameters in the theory. In this 
paper, these concepts have been applied to the effective potential, in application to inflation and 
quintessence models. This follows previous treatments of decoupling in effective potentials, ap- 
plied to electroweak physics [2lM23| . What is proposed here is that in interacting inflation or 
quintessence models, if there are disparate mass scales in the theory, then a mass independent 
renormalization scheme, such as the MS or MS schemes is not adequate in capturing the physics 
of decoupling. Rather the mass dependent renormalization scheme is more appropriate. In this 
scheme, as shown in this paper, heavy fields coupled to the inflaton or quintessence field have their 
quantum effects power suppressed in the effective potential. To implement the MDR scheme, the 
choice of renormalization scale, and thus the division between low and high energy scales in the 
effective potential, must be determined. 

Although the effective potential must be independent of renormalization scale, the issue is at 
what scale should the parameters in the theory be initialized. As noted in this paper, this is an 
underlying ambiguity in inflation and quintessence models and can lead to differing results. We 
discussed two options which we called the low-energy and high-energy approaches in Sect. IIVI 
In the low-energy approach the parameters are initialized at a renormalization scale at order the 
lowest mass scale in the system. One consequence of this approach is, if the inflaton or quintessence 
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field is the field with lowest mass, and all other helds interacting with it are much heavier, then 
there will be negligible quantum corrections, due to the decoupling effects. Alternative in the high- 
energy approach, the parameters are initialized at some high energy scale, possibly if the theory 
were embedded in some higher theory and some symmetry principle determined the parameters 
at high scale. In this case, one can use the renormalization group equations to evolve the scale to 
where one wants to do physics. Thus in the absence of the inflation or quintessence model being 
embedded in some higher theory, if the model has disparate mass scale, the low-energy approach 
is technically a viable option, in which case even if the model is moderately interacting with other 
fields, because of the decoupling effect it can still produce the very flat potentials needed. 

In cosmology there has been a common practice in the treatment of the scalar field effective 
potential, whereby the effect at one-loop order of any quantum field coupled to this scalar field 
leads to a Coleman- Weinberg correction term in the quantum corrected effective potential @. 
However this procedure is not always applicable. In particular, when the mass of a quantum field 
is greatly in excess to all other masses in the system, the decoupling theorem implies the quantum 
corrections from this heavy field are suppressed in powers of the light-to-heavy mass ratio. 

There are many common models in inflation and quintessence where these suppression effects 
become valid. For inflation, in all large field models, often terms chaotic inflation models [2|, as 
well as some hybrid models fl9j inflation occurs when the inflaton background field value is very 
large (0) ~ m p . Thus any scalar or fermion field coupled to the inflaton with at least moderate 
coupling will have a very large mass, much bigger than the inflaton mass. The quantum correction 
from such fields will thus be greatly suppressed during the inflation period due to decoupling. Also 
in small field models [2(| after inflation, as the inflaton background field value grows, it is possible 
fields coupled to the inflaton acquire large mass and thus there quantum corrections might become 
suppressed. In quintessence models, the dark en ergy regime in most models occurs when the scalar 
field background value is very large <j> ~ m p 0, fill [38j . Thus once again during the quintessence 
regime, quantum corrections from other scalar and fermion fields coupled to the quintessence field 
will be highly suppressed from decoupling effects. 

This decoupling effect can have significant importance to the building of inflation and 
quintessence models. For inflation it implies the inflaton field can be coupled more strongly to 
other fields. Stronger couplings can have beneficial effects in leading to more robust reheating 
after inflation or in the case of warm inflation models 38] , can lead to a greater parameter regime 
and many more possible models. Moreover inflation models in general can be less dependent on 
supersymmetry. In quintessence models, the potential during the dark energy regime typically 
has to be so flat, that the quintessence field often is just added as an additional field without any 
other dynamical purpose aside from driving the late time dark energy expansion of the universe. 
However the decoupling effect means that at early times the quintessence field might be interacting 
with other fields and produce dynamical consequences and at later time as the quintessence field 
background value becomes large, these other fields become very massive, thus decoupling, which 
then leads to a almost noninteracting field with the required ultra-flat potential to drive the late 
time dark energy phase. 
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Appendix A: Mass Dependent Renormalization Scheme 

In this Appendix we summarize the mass dependent renormalization scheme [24| used to obtain 
the RGEs for mass parameters and couplings Eqs. (PI51) - (|M)) . In order to set the notation and the 
procedure, we start by considering the simplest example, a single scalar field with potential: 



(53) 
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(a) [b) 

FIG. 6: (a) One-loop scalar self-energy diagram, (b) One-loop correction to the proper scalar quartic 
vertex. 

The renormalized Lagrangian including the counterterms is then: 



= ^d,m - -^mW - ^A^ 4 - Z^Q , (54) 



2 " T ' 2 f 4! 



(/>(] 



where 0, m 2 and A^ are (finite) renormalized parameters, and </>o = ^ TO 
A^o = Z^ 2 Zx^Xc/,, the (infinite) bare parameters, and 

6zt = Zt-l, (55) 

being the wave function renormalization constant. 

Loop corrections are computed with the Lagrangian given by the first line in Eq. (j5"4")l . includ- 
ing the corrections given by the counterterms. The first step is then to regularize the divergent 
integrals, and for that we use dimensional regularization: evaluating the integrals in d = 4 — e 
dimensions, and then taking the limit e going to zero. The divergent term at one-loop is isolated 
as the single pole when d — 4 (2/e term). The renormalized (finite) mass and coupling are defined 
by imposing suitable normalization conditions on the n-point 1PI Green functions r' n ) at some 
arbitrary scale p 33]. The relation between the bare and renormalized n-point 1PI functions is 
given by: 

r ( "V) \ P *=^ <M"V) , (56) 

where as before the subscript "0" denotes the bare quantity. The normalization condition fixes the 
counterterms that cancel out the divergent terms; this is equivalent to define the renormalization 
constants, given the relation between these and the counterterms introduced in section|TTl Eqs. (j2"2l) - 
We now derive them explicitly in the MDR scheme. 

The 2-point 1PI renormalized function F^ 2 \ including the contributions from the counterterms 
and the radiative correction II (p 2 ) (Fig. (JSJa)), is given by: 

r( 2 >( P 2 ) = P 2 - m 2 , + 5 z y - Sm% + ^ ] {p 2 ) ■ (57) 

The counterterms, or equivalently the renormalization constants, are fixed by demanding L^ 2 ' to 
be that of a free-field theory with running mass parameter m^(/i) [l8| . at the renormalization 
scale fi: 

T( 2 Hp 2 )\ p ,^(p 2 -m 2 (p 2 ))\ p2=fl2: (58) 

where p 2 is the incoming euclidean momentum 3 . The one-loop contribution is given in this case 
by: 

I#V) = ^L( m 2 / A 2 ) TO 2 , (59) 




L K//n = sh + i-^i , (60) 



In this Appendix p 2 will refer to the Euclidean momentum, and we drop for simplicity the subindex "E" hereon. 
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where the scale /t is introduced in the regularization procedure because of dimensional reasons, 
and 2/e = 2/e — ■fE + hi47r. When the scalar field does not couple to fermions, the wave function 
renormalization constant at 1-loop does not receive any contribution and therefore: 

Z^ = l. (61) 

On the other hand, the normalization condition fixes the mass counterterm: 

6m 2 = nW V), (62) 
and from Eq. (1221) we can read the mass renormalization constant: 

Z mr l-^L(m^ 2 ). (63) 

The relation between the renormalized mass parameter and the bare parameter is given by: 

ml{p) = Z^Z m yn\ , (64) 

and thus, taking the derivative with respect to the renormalization scale, one obtains the RGE for 
the mass parameter: 

* a In \i v a In fx 

Thus, given that at one-loop the radiative correction Ii^\p 2 ) is independent of the external mo- 
mentum, this (quadratically) divergent contribution can be reabsorbed into a redefinition of the 
mass parameter. 

For the renormalized vacuum energy, the situation is quite similar. The zero-point 1PI function 
at 1-loop are given by vacuum diagrams which do not depend on any scale except that of the mass 
of the particle in the loop. Therefore, 

4 

Z n Sl = SI - ^L(mJ/A 2 ) , (66) 

and 

Pa = . (67) 

The vacuum contribution will be fixed by the boundary conditions, say Sl(fi) — 0. Either at higher 
orders, or when the scalar couple to fermions, we will have Za, ^ 1 and (3 m 2 ^ 0, but still the pure 



vacuum contributions give /3q = 0. 

The RGE for the quartic coupling is derived in a similar manner, obtaining the coupling 
renormalization constant by imposing the normalization condition on the 4-point 1PI function. 
The renormalized 4-point 1PI function at 1-loop is given by: 

r( 4 V) = -A -5A + ^r( P W,) 

where T(p 2 ,m 2 ) is the contribution from the 1-loop diagram Fig. (|6lb): 



TV 2 2x 1 [2 f l . . , m 2 + x{l - x)p 2 l 

r(p ' m)= i6^u-y dxH — J- — } ) ■ (69) 



We impose the normalization condition: 

r^(p 2 )\ p2 ^-x^), (70) 

which defines the coupling counterterm: 

8\^^\lr(n 2 ,m 2 ), (71) 
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and using Eq. (|23[) . the coupling renormalization constant Z x : 



2 

The relation between renormalized and bare coupling is given by: 

72 r?-l 



\4,<jjL)=Z*Z^\ 40t (73) 
and as before, taking the derivative the RGE for the coupling reads: 

dX ^ - J^Li^aJ (74) 
dln M " (4tt) 2 21 0j ' 1 j 

= 1- ln^^5±l, (75) 

where = m^//i 2 . Threshold effects are included in the effective coupling through the momentum 
dependence p 2 = /i 2 of the radiative corrections 4 , such that the coefficients of the RGEs are 
modulated by a threshold function i*2(a). The latter reduces to one in the massless limit, a = 0, 
but it goes to zero when a = m 2 / goes to infinity. That is, in the massless limit one recover the 
same RGEs for the effective couplings than those computed in a mass independent scheme, like 
the MS scheme. But in the opposite limit, for a heavy state with o»l, the contribution is more 
and more suppressed as the ratio a increases. 

Having set the scheme in the simplest model, we can extend it now to the case of study, adding 
the second scalar field \i with potential: 



n x > = < i + \m\<? + ^ + im 2 x 2 + ^x 4 + \g 2 x 2 4> 2 ■ < 7(i ) 



Adding the corresponding counterterms, the renormalized Lagrangian for the scalar sector is then: 



+ % W + %e^x - ^ _ ^ _ ^ _ ^ _ ^ V _ sn . (77) 

We also consider Yukawa interactions between x an d Np fermions with the fermionic La- 
grangian given by: 

C F = *„(i7 M a M -m/)* a -^X*a* + ^/*a*-*W/* a *-J/lX*a*, (78) 

where we have taken for simplicity a common mass m/ and Yukawa coupling h for the fermions. 

The interaction term given by the coupling g 2 only adds a constant (quadratically) divergent 
term to 11^, through the same diagram than in Fig. (||)Ja) but with a x internal line. Therefore, 
still we have Z^ = 1 and 

m l Z ml = ™^- ^m\L(m\l ' p 2 ) ~ g 2 m 2 x L{m 2 J p 2 ) , (79) 

and thus /3 m 2 = 0. On the other hand, Z x and Z m 2 receive a p-dependent contribution from the 

loop of fermions. The diagrams contributing to the x ne ld 2-point function are given in Fig. J7]), 
with 

n W (p 2 ) = 2h 2 N F T( P 2 ,m 2 ) P 2 +8h 2 N F r( P 2 , m 2 f ) m 2 + ^-L{m 2 J fi 2 )m 2 x + g 2 L(m 2 / fi 2 )m 



4 In practice, the running effective coupling can be obtained by taking the derivative of the one-loop IPI Green 
functions with respect to the momentum, and then replacing p 2 = fj, 2 . 
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FIG. 7: One-loop self-energy diagram contribution to the \ propagator. Dashed lines represents scalars, 
and fermions are given by solid lines. 



M ■ B. 



(O) 



<t>,x 



X 



X s 



m 

FIG. 8: (a) One-loop proper vertex correction to A x . (b) One-loop proper vertex correction to g 



For the couplings, the renormalization constants Z x and Z m 2 are given then by: 

Z x = l-2h 2 N F T(^ 2 ,m 2 f ) , (81) 

2 



[Z m 2 = m 2 x - 8h 2 N F m 2 T(^ 2 , m 2 ) - ^-m 2 x L(m 2 x / (i 2 ) - g 2 m 2 L(m 2 / ft 2 ) . (82) 



The renormalization constant Z x ^ picks up a new term similar to that in Eq. (|72[) due to the 
loop with 2 massive x states, i.e., a diagram similar to that in Fig. (|B]b) but with x running in 
the loop. The 1-loop diagrams contributing to the X x and g 2 coupling renormalization constants 
are given in Figs. (|5Ja) and (j5]b) respectively. The set of renormalization constant are then: 

\Zll = A - h 2 r(v 2 ,m 2 ) - G^r^m 2 ) , (83) 
\ X Z^ = X x -^X 2 x T^ 2 ,m 2 x )-6g 4 T( f i 2 ,m 2 )-24h 4 N F r(p 2 ,m 2 ), (84) 

*V = .g^i-^r^^^-^r^^^-Vr,^ 2 ,^^^)) , (85) 



where: 



1 , m\x + iti 2 (1 — x) + x(l — x)p 2 
16ir 2 \e Jo A 2 



r 2 {p,m^mf) = T ^[--l dxln( - )). (86) 



To obtain the RGE for the Yukawa coupling we also need to renormalize the inverse fermion 
propagator S* -1 and the Yukawa vertex T^ 3 ' . The corresponding 1-loop diagrams are given in Figs, 
(jnia) and (jHJb). The normalization condition for S^ 1 and r' 3 ) are similar to Eqs. (|5"5|) and Eqs. 

S-W^ = (^-m f (p 2 ))\ p ^, (87) 
r( 3 W=M 2 = M/*). (88) 
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(a) (6) 
FIG. 9: (a) One-loop fermion self-energy, (b) One-loop Yukawa vertex correction. 

and from the above equations we can obtain Zf, Z^, and Z mj for massive fermions: 



2 



Z S = l-^T f (^,m^mj), (89) 
Z- 1 = l-h 2 r 2 (n 2 ,m 2 x ,m}), (90) 
Z m/ = l-h 2 r f (n 2 ,ml,m 2 f ), (91) 



where: 



i 2 2 2n 1 f 2 f 1 j , ,m 2 x + toKI — a;) + x(l — x)p 2 \ 
r^p 2 ,™ 2 ,™ 2 ) = dxa;ln(^ ^ ^ i '—)\ . (92) 

Finally, from the above renormalization constant, and the relations: 

(93) 
(94) 
(95) 
(96) 
(97) 
(98) 

the RGEs (|28l) - (f3"3"|) including the threshold functions are easily derived (setting m/ =0). The 
latter are given by: 

. . ,, . , dr 2 (p 2 , m 2 , 0) , 1 + a ,„ , 

Fx a = -{An) 2 2V ; ' >- = l-a\n , 99 

alnp^ a 

F 2 (a) = -^f^^ = 1 - -^=ln gg±l , (100) 







m 2 x (fi) 


= ^x Z m^ m x0 > 


Mm) 


= ^I^M: 




= Z x Z \ x ^xo> 


<?V) 


= Z^Z x Z g ^g 2 , 




= z x z f Z h' h o 



2 fdT f (p 2 ,m 2 ,0) , dr 2 (p 2 ,m 2 ,0)\ 1 , . . , 1 



F 3 (a) = -(4tt) 2 2 ' ' + 2 ^ ' 2 ' = 1 + 2a(l - (1 + a) In — ) , (101) 

\ amp' amp z J a 

where in computing F\(a) from the loop with one light and a heavy m x 3> we have set 
= 0. F 2 (a) is the threshold function for a scalar loop with 2 equal massive states, Fi(a) that of 
a scalar loop with one massless and one massive scalar state, and F 3 (a) that with massless fermions 
and one massive scalar. 



Appendix B: power suppression of two-loops coefficients 



At two-loop order, the RG-improved effective potential is given by the one-loop effective 
potential with running parameters evaluated using the two-loop RG equations. Following the 
MDR prescription, again the optimal choice to fix the renormalization scale is below all massive 
thresholds; thus the two-loop effective potential reduces to the tree-level potential plus the two-loop 
RGE functions. Decoupling will be included in the latter through threshold functions, similarly to 
the one-loop RGEs. 
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At two-loop order, both the wave function renormalization constant Z$ and the mass param- 
eter one Z m i get a \i dependent contribution from the sunset diagram in Fig. (jlOla). whereas the 
quartic coupling A^ correction comes from Fig. (fTUl b). The two- loop beta functions in a mass 
independent scheme can be found for example in Refs. [30l . [joj . where one only would need to 
extract the divergent contributions from those diagrams. In the MDR scheme we need to carry out 
the full calculation of the diagram keeping the finite contributions. We will not attempt such a full 
two-loop calculation here, and we only want to argue that such diagrams gives a power suppression 
of the corresponding RGE coefficients that go at least like 0(/j, 2 /M 2 ) when M 2 ^> /i 2 , M 2 being 
the heavy mass running in the loop. For the first diagram in the vertex correction in Fig. dlOl b). 
this can be viewed as the 1-loop vertex correction but with the LHS interaction replaced by an 
effective (momentum dependent) vertex again like that of Fig. (JSlb). According to the one- loop 
calculation, this will give the corresponding suppression by the mass running in the loop also in 
the two-loop coefficients. The contributions of the order of 0(\,pg 4 ), and 0(g 6 ), will be therefore 
suppressed by the heavy mass M 2 . That of the order of 0(A?) is suppressed by factors (9(^ 2 /Af|) 
when /i <C M<j,. In the cosmology models studied in the text of this paper, generally ~ /i and 
so these diagrams are not suppressed due to decoupling effects, but rather because is always 
tiny. The second diagram in Fig. (fTUl b) gives a term (0(h 2 g A ) coming from the insertion of the 
fermion loop in one of the internal x lines. Again, this contribution will be suppressed by a factor 
0(/i 2 /M 2 , similarly to the 1-loop vertex correction without the fermion insertion. 

We will therefore concentrate on the sunset diagram in Fig. (jlOl a). and in particular on the 
contribution to the wave function renormalization constant. In general for the sunset diagram with 
three different internal scalar masses we have: 

nt(p 2 ,m 2 ,m 2 ,m 2 ) = ~ I V ,,,7 . i, „r. ,nj ) + /r/V- n, f) I . (102! 




(4^)%V,m 2 ) - -2{- + -(--\n^- lE + I m {p\m 2 )) 1 (103) 

(An) 4 B(p 2 ,m 2 ) = h- - lE + \ - 2 f dx [ dy(l - y) ln(M 2 (x, y) + p 2 y(l - y))) , (104) 
1 e 1 Jo Jo 

where A/ 2 (x, y) = (m\x + m|(l — x))y/ (x(l — x)) + m|(l — y), and A-j is a general quartic coupling 
with the appropriate symmetry factors. In particular, when we have three <fi running in the loop 
then Xf — A^, and with two x an d one <j> then A 2 = 12g 4 . The function I m (p 2 ,m 2 ) is the finite 
contribution to the mass renormalization. By adding the sunset contribution to the 2-point function 
Eq. (|57|) . the normalization condition Eq. (|58p fixes the non- vanishing two- loop wave function 



(2) (2) (2) 

counterterm S z , which defines the wave function renormalization constant Z^ ' = 1 + S z , and 

then the anomalous dimension of the field: 

(2) d\nZ { * ] 1 A 2 

= d\nn 2 = (4^j4^T2 F22 ( a<# " a *' a *^ +9 F ^( a x, a X' a 4>)) , (105) 

with the threshold function given by the parametric integral: 

^22(01,02,03) = 2 / dx f dy V ^~ V }* r, (106) 

Jo Jo o(x, y) + y(l — y) 

where a(x, y) = M 2 (x, y)j /1 2 . From the above expression, one can check that indeed when we have 
only massless states running in the loop then ^22(0) = 1, but when any of the masses a; goes 
to infinity, then ^22(04) goes to zero as 0(l/<Xi). The integral can be evaluated numerically for 
arbitrary mass parameters, and it behaves similarly to the 1-loop threshold functions in Fig. (JXJ) . 

From the renormalization condition Eq. (|95p . the two- loop beta function for A^ gets contri- 

(2*) (2) (2) 

butions from the anomalous dimension 7^ and the proper vertex = — d\nZ x ^ /dln^j,, 
4? = 4A 07 f> +7 ( 2 > 



1 



(47r) 4 \^(3^ 22 ( a *' a <t>} ~ 6( ^( a 0' H)) 
+A 05 4 (4F 22 (a x , a x , a ) - 24G(a , a x ) - 96 5 6 G(a x , a x ) - 432N F h 2 g 4 G F {a f , a x )) .(107) 
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FIG. 10: (a) Two-loops scalar wave function renormalization diagram. (b) Two- loops correction to the 
proper scalar quartic vertex. 



We have not computed explicitly the 2-loop vertex threshold functions G(cn, aj), Gp(ai, ah), but as 
argued above, we can expect them to have the correct limits when a* ^> 1, and therefore decoupling 
in the sense of power suppression is also maintained in the MDR scheme at the 2-loop level. 



Appendix C: 1-loop effective potential within the mass dependent renormalization (MDR) 

scheme 



The 1-loop radiative correction to the effective potential, computed using dimensional regu- 
larization, is given by: 

AV ^ = ( M *( L W/* 2 ) + \) + KMK/P) + 5)) ■ ( 108 ) 

The renormalized 1-loop effective potential is obtained by adding to AVreg the appropriate coun- 
terterms: 

AVW = AyW + ( Z ^ + \ml(Z m l - 1)0 2 + ^-(Z^ - 1)0 4 , (109) 

and plugging = 1 and the renormalization constants given in Eqs. (|79p. (|83|) . computed in the 
MDR scheme, one has 5 : 

» m if* M l 9 9 M l a M l 9 , 9 M l 

+ ^(ln-^-/K/ Ai 2 ))+.gV 4 (ln-A-/K/ M 2 )) , (110) 

where 



VI + 4a- 1 



1(a) = dxln(a + x(l - x)) = In a - 2 - yl + 4a In -^=== . (Ill) 

Jo V 1 + 4a + 1 

Because the one-loop renormalization constant have been computed in the previous section in the 
symmetric phase of the theory (i.e., taking the vev of the field to vanish), the threshold functions 
1(a) depend on the mass parameters m 2 , instead of the physical effective masses M 2 relevant for 
the effective potential. Then, this expression would only lead to decoupling of heavy states when 
rritp ,m x ^> 4>, but not when tfi ^> m x . For example by taking \jl <C m<f>, m x , and expanding 
the threshold function 1(a) when a 1, we have: 

+^r(-^+-)+^(-^ + ■■•)), (112) 

4 ora^ om^ y 



5 The cosmological constant is renormalized by demanding V(<f> = 0) = 0. 
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and unless ~ m x , we are still left with potentially large logs, ln</> 2 /m 2 , hK/> 2 /m 2 , and the 
original problem addressed in Refs. [U [25j|. On the other hand, in the particular limit that the 
mass parameters vanish, — m x = 0, one just recover the expression for the effective potential 
computed in a mass independent scheme: 



In this case, given that both mass scales are set by the vev of the field (f>, large logs could be 
controlled and resummed by taking for example fj, ~ cf>. 

This apparent failure of the MDR scheme can be related to the fact that the effective poten- 
tial is computed in the non-symmetric phase of the theory, after shifting the field <fi by its vev. 
One should therefore also impose the renormalization conditions and get the counterterms in this 
phase. After shifting the field, the propagators running in the loops depend now on the effective 
mass and, by repeating the calculation done in the previous section, one can derive similarly the 
renormalization constants now with the threshold functions depending on M a . Therefore, this 
accounts to replace m 2 by M 2 in both the logs and threshold functions in Eq. (|110[) , and the 
effective potential is then given by: 

^ = ^(ln^-/(M> 2 )) + ^(ln^-/(M^ 2 ))J . (114) 

This is similar to the effective potential obtained in Ref . |23| ■ In that work, decoupling is introduced 
in the effective potential through step-functions at each physical threshold Mj = fj,, Eq. (|20|) . In 
our case it is implemented through the threshold function I(at) obtained when computing the 
I-loop radiative corrections. 

Finally, the RG-improved effective potential is given by absorbing the log dependence on the 
running parameters, such that the I-loop potential is given by the tree-level potential evaluated at 
the boundary t* = In /!*//!: 

V eff = \ml{U)4> 2 {U) + ^^{U). (115) 

This can be shown by explicit integration of the RGEs, plugging the result in the tree-level poten- 
tial. At 1-loop order, the mass parameter and the field do not run, m^(t r ) — m^, and </>(i*) = 4>, 
and integrating the RGE for A^ in the MDR scheme, Eq. (|2"5)) we have: 

A^) * A^( M ) + (l{Ml/£) - I(M|/m 2 ) - In ^\ 

( J (K/it) - ^ M V^ - ln f ) ■ ( 116 ) 

Taking ^* -C M^, M x , this gives: 



- A,( M ) + ^ (in ^ - KMl/A + ^ (in ^ - 7(M 2 / M 2 ) ) , (117) 



and therefore: 




M M 
(In -J- - I(M> 2 )) + , 9 V(ln -f - I(M*/^)) ] , (118) 
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where = A^(/i) and g = g (/i). 
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